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Transverse momentum spectra of identified particles in high energy 
collisions with statistical hadronisation model 
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A detailed analysis is performed of transverse momentum spectra of several identified hadrons in high energy 
collisions within the framework of the statistical model of hadronisation. The effect of the decay chain following 
hadron generation is accurately taken into account. The considered centre-of-mass energies range from ~ 10 to 
30 GeV in hadronic collisions (7rp, pp and Kp) and from ~ 15 to 45 GeV in e + e~ collisions. A clear consistency 
is found between the temperature parameter extracted from the present analysis and that obtained from fits 
to average hadron multiplicities in the same collision systems. This finding indicates that in the hadronisation, 
the production of different particle species and their momentum spectra are two closely related phenomenons 
governed by one parameter. 



1. Introduction 

The idea of a statistical approach to hadron 
production in high energy collisions dates back 
to '50s 0] and '60s and it has been recently 
revived by the observation that hadron multi- 
plicities in e + e~ collisions agree very well with 
a thermodynamical-like ansatz [^JJ] . This finding 
has also been confirmed in hadronic collisions and 
it has been interpreted in terms of a pure statis- 
tical filling of multi-hadronic phase space of as- 
sumed pre-hadronic clusters (or fireballs) formed 
in high-energy collisions, at a critical value of 
energy density ||. In this framework, tempera- 
ture and other thermodynamical quantities have 
a statistical meaning which does not entail the 
existence of a thermalised hadronic system on an 
event-by-event basis. Stated otherwise, statistical 
equilibrium shows up only when comparing many 
different events, whilst in each of them the Gibbs 
law of equally likely multi-hadronic states applies. 
So far, this proposed statistical cluster hadronisa- 
tion model has been mainly tested against mea- 
sured abundances of different hadron species for 
a twofold reason. Firstly, unlike momentum spec- 
tra, they are quantities which are not affected by 
hard (perturbative) QCD dynamical effects but 
are only determined by the hadronisation process; 
indeed, in the framework of a multi-cluster model, 



they are Lorentz-invariant quantities which are 
independent of the cluster's overall momentum. 
Secondly, they are quite easy to calculate and pro- 
vide a very sensitive test of the model yielding 
an accurate determination of the temperature. 
However, in order to establish the validity of this 
approach, it is necessary to test further observ- 
ables and to assess their consistency with the re- 
sults obtained for multiplicities. Certainly, one of 
the best suited observables is the transverse mo- 
mentum of identified hadrons, where transverse is 
meant to be with respect to beam line in high en- 
ergy hadronic collision, and thrust or event axis 
in high energy e + e~ collisions. Indeed, such pro- 
jection of particle momentum is supposed to be 
the most sensitive to hadronisation or, conversely, 
the least sensitive to perturbative QCD dynam- 
ics. 

Actually, it has been known for a long time 
that transverse momentum distributions are 
Boltzmann-like in hadronic collisions and this 
very observation was pointed out by Hagedorn 
as a major indication in favour of his statisti- 
cal model of hadron production 0. It must 
be emphasized that the prediction of a thermal- 
like shape in principle only applies to particles 
directly emitted from the hadronising source, 
whereas measured spectra also include particles 
produced by decays of heavier hadrons. However, 
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most analyses do not take into account the dis- 
tortion of primordial hadronisation spectrum due 
the to hadronic decay chain and try to fit the 
data straight through it. This problem has been 
discussed in literature @ and an analytical calcu- 
lation has been developed to take into account the 
effect of two and three body decays |||| , which 
has then been used both for pp ||] and heavy 
ion collision |£-10 including most abundant res- 
onances. In this paper we introduce a method al- 
lowing to rigorously and exhaustively determine 
the contribution of all particle decays. Hence, by 
taking advantage of this technique, we have per- 
formed an analysis of many measured transverse 
momentum spectra of identified hadrons in a wide 
range of centre-of-mass energies for several kind 
of collisions. 

2. Statistical hadronisation and transverse 
momentum spectra 

The statistical hadronisation model |l2| as- 
sumes that in high energy collisions, as a con- 
sequence of strong interaction dynamics, a set of 
colourless clusters (or fireballs) is formed having 
certain values of mass, volume, internal quantum 
numbers and momentum, the latter being inher- 
ited from the hard stage of the process. Those 
clusters are assumed to give rise to hadrons ac- 
cording to a pure statistical law in the multi- 
hadronic phase space defined by their mass, vol- 
ume and quantum numbers. This approach dif- 
fers from another popular cluster hadronisation 
model [|l3| mainly because it gives clusters a vol- 
ume so that hadron production is ruled by the 
properly understood phase space rather than rela- 
tivistic momentum space. In this framework, the 
use of statistical mechanics and thermodynamical 
quantities, such as temperature, which need spa- 
cial dimension besides momentum space in order 
to be meaningful, is allowed. We emphasize once 
more that the introduction of such quantities does 
not entail any thermalisation process of hadrons 
after their formation, nor the existence of a ther- 
malised system event by event. 
Although many clusters with different momenta, 
volumes, masses and quantum numbers are 
formed, it can be shown that the average val- 



ues of many observables, e. g. particle multi- 
plicites, are the same as those relevant to one 
equivalent cluster having suitable values of vol- 
ume (namely the sum of all cluster volumes mea- 
sured in the rest frame of the equivalent cluster 
itself), mass and quantum numbers (namely the 
sum of all cluster quantum numbers) . The proof 
of this statement a lengthy one, requires the 
assumption of special probabilities governing the 
fluctuations of cluster masses and quantum num- 
bers for a given set of volumes. If the volume (or 
mass) of the equivalent cluster is large enough, it 
is then reasonable to take a canonical approach 
(i.e. introducing a temperature) in order to calcu- 
late mean quantities, instead of carrying out an 
involved microcanonical calculation. Therefore, 
even though actual clusters are small sized and 
microcanonical calculations would be needed to 
determine mean quantities within each of them, 
the choice of suitable mass fluctuation probabili- 
ties for the clusters allows one to calculate overall 
means dealing with only one large global struc- 
ture and with much fewer parameters. Arguing 
the other way around, it is not difficult to be 
convinced that the fluctuation probabilities to be 
chosen in order to achieve such reduction of the 
problem are exactly the same as those of obtain- 
ing a given set of masses by splitting a large clus- 
ter into a given number of clusters with volumes 
(Vi, . . . , Vn). This statement extends a reduction 
procedure to the microcanonical case which was 
proved in the canonical case [||, where, from the 
very beginning, clusters were given temperature 
and volume instead of mass and volume. Never- 
theless, it is not obvious how large the equivalent 
cluster should be for the canonical approxima- 
tion to apply. For the present, we have adopted a 
simple-minded a posteriori method consisting in 
justifying the canonical framework by its capa- 
bility of accordance with the data. A simple ar- 
gument to support the canonical approximation 
even at moderate values of volumes is the very 
large number of states (O(10 2 )), i.e. hadrons and 
resonances, which can be excited in a hadron gas. 
As far as single-particle transverse momentum 
spectra are concerned, a similar reduction theo- 
rem from many clusters to one equivalent cluster 
in the averaging procedure applies fl4|| , though 
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only approximately. In general, it can be shown 
that the primary spectrum of j th hadron species 
depends on transverse four-velocities UTi = 
flrjTi Q of the N clusters: 
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where /(mti, • ■ ■ , utn) is the transverse four- 
velocities distribution function. If one expands 
all single-cluster spectra dnj / dpT\i{uTi) in se- 
ries of UTi starting from a common value (ut) 
for all clusters, it can be proved that, at the 
zeroth order, the spectrum in eq. ([!]) becomes 
the same as that obtained for the aforementioned 
equivalent cluster, endowed with a tranverse four- 
velocity (ut)- This reduction possibly allows 
taking the canonical approach in order to calcu- 
late the transverse momentum spectrum since the 
equivalent cluster size is much larger than single 
cluster's: 
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where mx = 



Pt 



T is the temperature, 



and Ki , I are modified Bessel functions. Eq. (§) 
is the Boltzmann limit of quantum statistics and 
it is a very good approximation for all hadrons 
except pions [||. For the zeroth order approx- 
imation in the ut% expansion to be sufficiently 
accurate, the involved transverse four-velocities 
should be small, i.e. /(uti> • • • , Utn) — * al- 
ready for un "C 1. Since it is possible to 
choose the starting point (ut) of the series ex- 
pansion to make the first order term vanishing 
PI , the largest neglected term turns out to be 
0((uTi — (ut)) 2 ) which is very small provided 

1 In fact, we mean by transverse four- velocity the module 
of the spacial part of a velocity four-vector u = (7, f3y) 
with vanishing longitudinal component. 



that the aforementioned condition on / is met. 
The analysis presented in this paper assumes the 
validity of this approximation and primary spec- 
tra have been calculated according to eq. (||) . 
Even though eq. @ looks like the spectrum from 
a thermalised source with a superimposed flow, 
a popular formula in the heavy-ion community, 
it is well worth emphasizing that in fact this for- 
mula has nothing to do with flow. Indeed, the av- 
erage transverse four-velocity (ut) is not meant 
to be the mean value of cluster velocity distri- 
bution in a single event; rather, it is the aver- 
age transverse four-velocity of all clusters over all 
collision events. In other words, (ut) may well 
arise from single clusters emitted at high px in 
some events, e.g. following a hard parton scatter- 
ing, and does not imply by any means a collective 
event-by-event expansion. 




Figure 1. Measured spectra of six identified par- 
ticles in pp collisions at \/s = 27.4 GeV |2f|. The 
measured dn/dp^ spectra (above) have been nor- 
malised so as to have the same value at pr = 0. 
While slopes of dn/dp 2 ^ spectra are different, 
those of dn/drriT (below) are approximately the 
same for all identified particles. 
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3. Method of data analysis 

A peculiar prediction of a statistical picture in 
high energy collisions which is relevant to trans- 
verse momentum spectra is the so-called mj- scal- 
ing: there should be an apparent common slope 
for dn/dniT spectra of identified hadrons at a 
given centre-of-mass energy. By looking at eq. (Q) 
one can easily realize that this holds in the limit 
(ut) — ► 0. In principle, tut scaling applies only 
to primary hadrons, namely those directly emit- 
ted from the hadronising source. On the other 
hand, as stated in Sect. 1, most observed or recon- 
structed particles in experiments arise from de- 
cays of heavier hadrons. Those secondary decays 
may well distort the primary spectrum shape, 
thus spoling mj scaling. Nevertheless, it can be 
seen in fig. (|]) that bit scaling apparently also 
holds for measured final hadrons, at least in the 
examined collision, implying seemingly little dis- 
tortion from primary to final spectrum. Settling 
this and other related issues demands a thor- 
ough analysis taking into account the effect of 
decays. Therefore, measured pt spectrum of the 
j th hadron species should be compared with the 
sum of its primary spectrum and the contribution 
arising from all heavier hadrons decaying into it: 
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Whilst dnj/dp T | primary is given by eq. (§), the 
calculation of the second term in the right hand 
side of eq. (^) within a statistical-thermal frame- 
work is a very complicated problem as the num- 
ber of contributing resonances and decay modes 
is huge and also owing the presence of cascade 
decays (e.g. rj' — > p — ► 7r). In rcfs. p|-|lO[] an ana- 
lytical approach has been taken including two and 
three body decays of directly emitted resonances. 
On the other hand, we have adopted |HJ a mixed 
numerical-analytical method which allows a very 
accurate and exhaustive calculation. Our method 
can be summarised as follows: for any couple of 
hadrons (J, k), where k decays into j + X, either 
directly or through a cascade process, we have 
determined the momentum spectrum dnj/dp\ ^ 
in the equivalent cluster's rest frame via a Monte- 
Carlo method by randomly generating 200,000 k's 



decays, according to phase space, by picking k's 
three-momenta from a thermal distribution with 
a given temperature T: 
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Then, we have boosted the obtained drij / dp\p_Q 
distributions in order to compute the 
dn,-/d|>| J for a given transverse four- 
velocity (it-r) of the equivalent cluster, to be 
plugged into eq. (||), with an integral formula 
derived by the authors Q : 
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where F is the elliptic integral of the first kind 
and: 
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The main advantage of formula (||) is to allow a 
quick computation of transverse momentum spec- 
tra associated with particle decays for any tran- 
verse four-velocity once the corresponding mo- 
mentum spectra at zero velocity are known. To 
obtain it, we took advantage of the isotropy of 
dnj/dp|^Q distributions, which holds as long as 
decay products distribution is isotropic, which is 
true for simple phase-space decays with no po- 
larisation. The integration in the variable p' in 
eq. (^) has been performed numerically. The 
Monte-Carlo generation of dn^/dpLZo spectra, 
depending on the temperature, has been per- 
formed for all unstable hadrons and resonances 
up to a mass of 1.8 GeV and it has been repeated 
for temperature values between 140 and 190 MeV 
in steps of 1 MeV, so as to have a large enough 
range of temperatures where to search for minima 



5 



Table 1 

Results of the fit to hadron average multiplicities for the chosen set of hadronic and e + e~ collisions. The 
X 2 value in pp collisions is very high due to lack of systematic errors fl6|| in the measurements of particle 
multiplicities 



collision 


v^(GeV) 


T (MeV) 


VT 3 


(SB) 


X 2 /dof 


K+p 


11.5 


176.9±2.6 


5.85±0.39 


0.347±0.020 


68.0/14 


7T + p 


21.7 


170.5±5.2 


10.8±1.2 


0.734±0.049 


39.7/7 


K+p 


21.7 


175.8±5.6 


8.48±1.05 


0.578±0.056 


38.0/9 


PP 


27.4 


162.6±1.6 


14.17±0.66 


0.644±0.018 


313.9/29 


collision 


Vi(GeV) 


T (MeV) 


VT 3 


7s 


X 2 /dof 


e + e~ 


14 


167.3±10.4 


9.7±2.6 


0.795±0.089 


1.5/3 


e + e~ 


22 


172.5±5.3 


10.6±1.8 


0.767±0.091 


1.0/3 


e + e~ 


29 


159.0±1.8 


17.3±1.1 


0.710±0.039 


29.3/12 


e + e~ 


43 


162.5±5.4 


16.2±2.1 


0.768±0.064 


3.0/3 



in a fitting procedure. For any other temperature 
value, the spectra have been calculated by means 
of a linear interpolation. Overall, for hadronic 
collisions, we have generated about 1200 spec- 
tra associated to 144 unstable hadrons for each 
T value. 

The abundances of all hadrons, which are needed 
to assess the contribution of secondaries in the 
spectra of eq. (^), have been preliminarly set to 
the model values calculated by fitting the pa- 
rameters T, V and (ss) or 75 to the experimen- 
tally measured yields. This procedure amounts 
to exclude the overall normalisation dependence 
of drij /&pT\ k ~*3 on T in the fit so to keep only its 
shape dependence on T. 

Whilst the 75 parameter has already extensively 
been used in multiplicity fits , a new param- 
eter (ss) is introduced here for the strangeness 
suppression. The particle yields are computed 
by constraining the number of newly created 
ss pairs to fixed integers which fluctuate accord- 
ing to a Poissonian distribution. The mean value 
of the Poissonian distribution, namely the mean 
value of newly created ss pairs, is the adjustable 
fit parameter (ss) . For the present, this new 
parametrisation has been used only for hadronic 
collisions (see Table 1) because e + e~ collisions 
present difficulties related to the production of 
heavy flavoured quarks which makes calculations 
extremely slow. The mathematical features of 



the (ss) parametrisation are described in detail 
in ref. jTi) . 

The fitted values, which update those in refs. j|,[3) 
for hadronic and e + e _ collisions are quoted in Ta- 
ble 1 The fitted temperature in pp collision in 
the present analysis is somewhat lower compared 
with ref. [ p| o wing to an updated set of hadron pa- 
rameters [ |19| , an upgraded procedure of x 2 min- 
imisation also taking into account correlations in 
the effective variance method pOfl , as well as the 
replacement of 75 with (ss) . The fitted tempera- 
ture values show a very interesting trend consist- 
ing in a slight but significant increase as centre- 
of-mass energy diminishes; this is in agreement 
with the observation pointed out in ref. pl[ . 

4. Data set and results for hadronic colli- 
sions 

The collisions and relevant centre-of-mass en- 
ergies to be analysed have been primarily cho- 
sen on the basis of a fairly large number of mea- 
sured particle yields and spectra. Particle multi- 
plicities are needed to establish with some accu- 
racy the contribution of secondary decays while 
many identified hadron spectra allow performing 
a test of universality on the slopes governed by 

2 The measurements used to perform this fit have been 
gathered through the Durham reaction data database ]f8| ; 
the relevant references will be quoted in ref. JL4{ . 
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Figure 2. \ 2 contour plot in the T — (ut) 
for 7r + p\ spectrum in pp collisions at yfs 
GeV. 



plane 
= 27.4 



mental errors and errors arising from the uncer- 
tainty on masses, widths and branching ratios 
of hadrons decaying into the examined hadron, 
according to the effective variance method (see 
ref. [0 for a detailed description). Once the 
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T and (ut)- This criterion led to singling out 
mainly four collision systems for hadronic colli- 
sions: K+p at 75 = 11.5 p§ and 21.7 GeV (23), 
7T+p at Jl = 21.7 GeV @ and pp y/s = 27.4 
GeV [^5|. For each identified particle transverse 
momentum spectrum, three free parameters have 
been determined by fitting eq. (J3J) to the data: 
T, (ut) and an overall normalisation parameter 
A. The fits have been performed in the variables 
Pt or p\ according to the available experimental 
spectrum; henceforth, we will denote by dn/dpr 
both dn/dpr and dn/dp^. As we have mentioned 
before, overall multiplicities of hadrons in eq. (||) 
have been fixed to model values by using fit pa- 
rameters in Table 1. 

The fit is performed by minimising the following 
X 2 : 



x 2 = E 



[ /bin, d P T d "/ d PT - (Exp.value), 



"(7) 



The errors Oi in the denominator of eq. (Q) have 
been taken as the sum in quadrature of experi- 



Figure 3. Contour plot in the T — (ut) plane 
of E l xl for 7T+p collisions ^ = 21.7 GeV. The 
identified particles involved in the sum are ir~ , 7r°, 
K^, p°, K*° and K*°. The two vertical lines define 
the 1-sigma band corresponding to the tempera- 
ture fitted by using average particle multiplicities. 
The Y^ii Xi quoted in the picture is its value in the 
local minimum pointed by the solid line. 



normalisation parameter A is fitted, a study of x 2 
contour lines is performed in the T — (ut) plane 
in order to look for possible local minima. In- 
deed, the presence of several local minima along 
a band is a quite general feature of the x 2 ( see 
fig. (||)) owing to the strong anticorrelation be- 
tween T and (ut)- Hence, different solutions are 
possible for the (T, (ut)) pair. In order to enforce 
universality of such parameters among different 
hadron species we have chosen a solution for each 
particle according to the following procedure: 
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Figure 4. Contour plot in the T — (ut) plane of 
*£,iXi for K+p collisions y/s = 21.7 GeV. The 
identified particles involved in the sum are 7r~, 
7T°, K%, p°, K*0 and f 2 (1270). The two vertical 
lines define the 1-sigma band corresponding to 
the temperature fitted by using average particle 
multiplicities. The J^. \i quoted in the picture 
is its value in the local minimum pointed by the 
solid line. 



Figure 5. Contour plot in the T — (ut) plane of 
J2i Xi f° r PP collisions y/s = 27 A GeV. The iden- 
tified particles involved in the sum are 7r + , 7r~, 
7T°, K+, K~, T], p°, and f 2 (1270). The two verti- 
cal lines define the 1-sigma band corresponding to 
the temperature fitted by using average particle 
multiplicities. The xf quoted in the picture 
is its value in the local minimum pointed by the 
solid line. 



1. firstly, we have determined the contour lines 
of YliLi Xi as a function of T and (ut) by 
fixing all normalisation parameters to their 
fitted values, where N is the number of dif- 
ferent hadron species spectra for a given col- 
lision. The minimisation of the sum of all 
\ 2, s amounts to fit a common value of T 
and {ut) for all measured hadrons. 

2. secondly, among the different local minima 
we have chosen the one lying in or closest 
to the band defined by T fitted with aver- 
age multiplicities along with its error (see 

ngs. dig). 

3. finally, all single particle spectrum fits have 
been repeated seeking a local minimum as 



close as possible to the previously deter- 
mined Xi minimum. 

It must be noted that this procedure aims at 
achieving the best agreement between tempera- 
tures determined by using different observables, 
namely particle yields and transverse momentum 
spectra. The 53«=i Xi contour plots are shown in 
figs. (H;§,|[) for the examined collisions. There is 
a good agreement between the T obtained from 
multiplicities and the location of minima in the 
T - {u T ) plane. For K+p and 7r+p at y/s = 21.7 
GeV the absolute minimum lies nearly about the 
centre of the band demonstrating an intriguing 
correlation between the slope of transverse mo- 
mentum distributions and the slope of average 
multiplicities as a function of mass. Moreover, 
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the (ut) values turn out to be small, thus jus- 
tifying the series expansion described in Sect. 2. 
Conversely, in K + p collisions at ^/s = 11.5 GeV 
a clear disagreement emerges for K + p at ^fs = 
11.5 GeV between differently determined temper- 
atures (see fig. (j|)) which can be explained with 
the inadequacy of our used parametrisation at low 
energy where exact pt conservation, neglected in 
the canonical framework, is expected to play a 
significant role. The summary plots of the sin- 
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Figure 7. Above: 7r + measured and fitted spec- 
trum in pp collisions at yfs = 27 .4 GeV. The lines 
lying below the top solid line represent the cumu- 
lative contributions to the fitted spectrum of sec- 
ondary 7r + 's arising from the decays of doubly and 
singly strange baryons, other baryons, strange, 
charged and neutral mesons. Below: residuals 
distribution. 



Figure 6. x 2 contour plot in the T — (ut) plane 
for 7T° and Kg p T spectra in K + p collisions at 
y/s =11.5 GeV. The two vertical lines define the 
1-sigma band corresponding to the temperature 
fitted by using average particle multiplicities. Lo- 
cal minima lie outside the 1-sigma temperature 
band. 

gle particle fits are shown in figs. (|[ ^,|l^) while 
an example of fitted spectrum in pp collisions is 
shown in fig. (0). In general, very good fits (low 
X 2 's) are obtained with good agreement between 
temperatures and average transverse boost veloc- 
ities fitted for different identified hadrons, espe- 
cially in pp collisions. Some discrepancies show 



up in K + p and 7r + p collisions where K*'s seem 
to have slopes steeper than expected. It must be 
noted that error estimates for T and (ut) are still 
very rough. 

5. Results for e + e collisions 

We have also performed a similar analysis for 
e + e~ collisions at centre-of-mass energies be- 
tween 14 and 44 GeV. Compared with hadronic 
collisions, the analysis of e + e~ data presents sev- 
eral additional difficulties. First of all, the pro- 
duction of heavy flavoured hadrons can not be 
neglected and the inclusion of many more decay 
channels is implied. A major issue is the fact 
that pr is defined with respect to the so-called 
thrust or event axis. This direction is not known 
a priori, unlike in hadronic collisions where pt is 
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Figure 8. Summary of the fits to transverse mo- 
mentum spectra of identified particles in pi + j> 
collisions at y/s — 21.7 GeV. From top to bot- 
tom: temperatures (the dashed line is the value 
obtained by fitting multiplicities), average trans- 
verse four-velocity (ut) and minimum x 2 /dof. 
The errors on T and (ut) are still very rough 
estimates. 



Figure 9. Summary of the fits to transverse mo- 
mentum spectra of identified particles in K + p 
collisions at yfs — 21.7 GeV. From top to bot- 
tom: temperatures (the dashed line is the value 
obtained by fitting multiplicities), average trans- 
verse four-velocity (ut) and minimum x 2 /dof. 
The errors on T and (ut) are still very rough 
estimates. 



defined with respect to the beam line, and must 
be determined on an event by event basis. Hence, 
whatever the algorithm being used, this very fact 
introduces a bias on transverse momentum spec- 
tra because momenta projection are used to de- 
termine the event or thrust axis itself; often, this 
is done by just minimising the sum of particle 
Pt's. In spite of this problem, we have repeated 
the same analysis performed for hadronic colli- 
sions by using measured charged tracks transverse 
momentum spectra |p6| with respect to thrust 
axis. The x 2 contour plots are shown in fig. (11) 
along with T bands obtained from multiplicity 
fits and resulting x 2 's minima. It can be seen 
that no minimum is found within the band for 
e + e~ collisions at \fs = 14 GeV, a situation look- 
ing very similar that found in K + p collision at 



y/s = 11.5 GeV. It is also found that fitted (ut) 
increases quite rapidly with centre-of-mass energy 
reflecting the rise of average pt of gluon radia- 
tion. Since our parametrisation of transverse mo- 
mentum spectra with an average transverse four- 
velocity (ut) requires (ut) to be <C 1 (see dis- 
cussion in Sect. 2), it is quite natural that fit 
quality deteriorates as centre-of-mass energy in- 
creases. This is indeed found already at 44 GeV 
where (ut) > 0.55 cannot longer be considered 
small. As energy increases, particle px spectra 
become sensitive to the shape of clusters ut dis- 
tribution (related in turn to radiated gluon pt 
spectrum) and not only to its mean value. Fur- 
ther on, they become dominated by the shape of 
gluon radiation spectrum whilst hadronisation px 
plays the role of a small superimposed noise. 



10 




Figure 10. Summary of the fits to transverse 
momentum spectra of identified particles in pp 
collisions at yfs — 27 A GeV. From top to bot- 
tom: temperatures (the dashed line is the value 
obtained by fitting multiplicities), average trans- 
verse four-velocity (ut) and minimum \ 2 /dof. 
The errors on T and (ut) are still very rough 
estimates. 



Figure 11. \ 2 contour plot in the T — (ut) 
plane for charged tracks pr or y\ spectra in 
e + e~ collisions at four different centre-of-mass 
energies ranging from 14 to 44 GeV. The two ver- 
tical lines define the 1-sigma band corresponding 
to the temperature fitted by using average par- 
ticle multiplicities. For e + e~ collisions at yfs = 
44 GeV, local minima lie outside the considered 
(ut) upper bound and the minimum x 2 is poor. 



6. Discussion and conclusions 

The analysis of transverse momentum spectra 
of several identified hadrons, accurately taking 
into account the effect of hadron decays, indi- 
cates that they can be well reproduced within the 
statistical hadronisation model for several high 
energy hadronic collisions. Furthermore, a good 
agreement is found between the temperatures es- 
timated by fits to average particles multiplicities 
and those extracted by fitting transverse momen- 
tum spectra. This can be seen especially by com- 
paring figs. (HJUJI); in fact, the location of the best 
T value for transverse momentum spectra, deter- 
mined by % 2 's lowest contour line, is apparently 
correlated to the best T value defined by parti- 
cle multiplicities as spotted by the vertical band. 
Indeed, those temperatures move toghether from 



~ 162 MeV in pp collisions at y/s = 27 A GeV to 
~ 175 MeV in 7r + p and K + p collisions at y/s = 
21.7 GeV. In our view, this finding is a strong in- 
dication in favour of one of the key predictions of 
the statistical hadronisation model, namely the 
existence of a close relationship between the laws 
governing the production of particles as a func- 
tion of their mass and, for each particle species, 
the production as a function of momentum (mea- 
sured in the rest frame of the cluster they belong 
to) at the hadronisation. 

Besides this main conclusion, there are many re- 
marks to be made. First of all, our calculations 
assumed a canonical framework and so they are 
expected to have only a limited validity at low 
centre-of-mass energies where the effect of exact 
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transverse momentum conservation must show 
up. Indeed, in e + e~ collisions at a centre-of- 
mass energy of 14 GeV and K + p collisions at 11.5 
GeV, a clear disagreement is found between the 
temperatures determined in the two fashions (see 
fig. (||)). Secondly, whilst the results obtained 
in pp collisions for different particles are in very 
good agreement among them, there are some sig- 
nificant discrepancies between different particles 
in 7r + p and K + p collisions, especially for K*'s, 
which are not understood at the present. Fi- 
nally, it is evident in the analysis of moderately 
high energy e + e~ collisions, at y/s = 44 GeV, 
that our parametrisation of transverse momen- 
tum spectra as a function of an average trans- 
verse four-velocity alone, is not accurate enough. 
As the average px of radiated gluons and, as a 
consequence, of hadronising clusters, rises, par- 
ticle spectra are influenced more by the shape 
of clusters transverse four- velocity spectrum than 
by primordial hadronisation px spectrum so that 
more accurate calculations are needed involving 
perturbative QCD to perform such an analysis. 
Finally, in the very high energy regime, px spec- 
tra become insensitive to hadronisation and, as a 
consequence, it is no longer interesting to study 
its properties. In summary, the presently used 
parametrisation of px spectra is found to work 
well only in a limited centre-of-mass energy range 
(roughly between 20 and 30 GeV). At lower en- 
ergies, the study of px is still of great inter- 
est in probing the statistical features of hadro- 
nisation but complex microcanonical calculations 
are required. At higher energy, hadronisation is 
swamped by hard QCD dynamics in px spectra 
and all it can do is to add a little smearing. 
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